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ABSTRACT 


The  detection  of  periodic  signals  of  arbitrary  wave  shape 
with  random  time  delay  in  additive  white  Gaussian  noise,  is  a 
problem  of  practical  significance  in  radar  and  communication 
applications . 

In  this  thesis,  the  analysis  and  design  of  optimum  and 
suboptimum  receivers  for  detecting  signals  as  described  above 
has  been  carried  out.  The  design  of  optimum  (in  minimum  proba¬ 
bility  of  error,  Pe  sense)  receivers  is  based  on  the  likelihood 
ratio  test  under  the  assumption  of  low  SNR  conditions.  The 
design  of  suboptimum  receivers  is  based  on  heuristic 
approaches  that  intuitively  yield  reasonably  good  performance. 
Examples  have  been  analyzed  in  order  to  present  numerical 
results  in  graphical  form  on  the  performance  of  the  receivers 
under  different  assumptions  of  wave  shapes  and  p.d.f.  on  the 
random  time  delay  associated  with  the  signal. 
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I.  INTRODUCTION 


Classical  noncoherent  signal  detection  is  generally  under¬ 
stood  to  mean  the  detection  of  a  sine  wave  with  random  phase 
or  time  delay  in  additive  white  Gaussian  noise  (WGN) .  However, 
the  extension  of  this  problem  to  a  general  noncoherent  problem 
involving  the  detection  of  ah  arbitrary  periodic  signal  with 
random  time  delay  in  additive  WGN,  has  received  little  attention 
A  case  of  particular  practical  interest  involves  the  detection 
of  a  baseband  square  wave  with  random  time  delay  in  additive  WGN 

In  digital  communication  systems,  the  message  to  be  trans¬ 
mitted  is  encoded  into  a  sequence  of  binary  digits.  Typically, 
these  digits  represented  by  the  logical  states  '1'  and  'O'  are 
transmitted  by  sending  a  suitably  chosen  set  of  pulses  which 
are  distorted  during  both  transmission  and  reception.  The 
effect  of  this  distortion  is  that  at  the  receiver,  it  is  no 
longer  possible  to  determine  exactly  which  waveform  was  actually 
transmitted.  We  can  model  the  transmission  as  well  as  the 
distortion  introduced  in  the  receiver  itself,  as  random  noise. 

In  all  the  analysis  carried  out,  it  will  be  assumed  that  the 
random  noise  can  be  modeled  as  additive,  white  Gaussian.  The 
problem  then  is  one  of  deciding,  on  the  basis  of  noisy  obser¬ 
vations,  whether  the  transmitted  waveform  corresponds  to  a 
logical  ' 1 '  or  a  logical  'O'. 

In  this  thesis,  the  analysis  and  design  of  noncoherent 
receivers  that  optimally  (and  in  some  cases  suboptimally ) 


It  is  clear  from  Equation  (2.30)  that  this  simple  receiver 
can  only  operate  provided  Vq  ?  0.  If  VQ  is  indeed  zero,  some 
other  approach  must  be  used  for  the  purpose  of  deciding  on 
which  hypothesis  is  true  upon  reception  of  r(t).  In  the  sequel 
we  discuss  other  detection  approaches  which  will  work  even  if 
V0  =  0. 

The  performance  of  this  receiver  is  evaluated  as  follows. 

Let 


JL 


r(t)  dt 


(2.31) 


Conditioning  on  A ,  we  see  that  £  is  a  conditional  Gaussian 
random  variable,  where 


eU|hq} 


0 


(2.32) 


EU|H1#  A} 


v ( t-A ) dt 


.2iTk,  T  .2-rTk. 

00  -]-m— A  O  }-S—t 

V  I  Vk  e  f  e  dt 

k=-oo  0 


A 

V  nT  =  m. 
0  o  1 


(2.33) 


Furthermore 


25 


However,  since 


so  that  the  test  can  be  expressed  as 


T 

o 

/  r ( t ) dt 

0 


> 


< 


(2.30) 


and  the  receiver  structure  takes  on  the  simple  form  shown  in 
Fig.  2.2. 


Figure  2.2  Nonchoherent  Receiver  (Uniform  p.d.f., 
2  Term  Approximation) 
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The  test  of  Equation  (2.23)  is  valid  provided  Equation 
(2.20)  holds. 

Observe  that  in  order  to  limit  the  approximation  error  to 
10%  in  Equation  (2.21),  |x|  £  0.39. 

In  the  remaining  analysis  we  assume  equal  prior  probabili¬ 
ties,  namely  P{HQ}  ~  \  -  P{H.}  and  cost  assignment  C. .  =  1  -  6. ., 
so  the  threshold  of  Equation  (2.23)  becomes 

N  E/N 

n*  =  -fie  °  -  1]  (2.24) 

In  order  to  further  analyze  the  test  of  Equation  (2.23),  we 
consider  the  following  two  cases: 

1.  Uniform  p.d.f.  on  X 

2.  Nonuniform  p.d.f.  on  X 
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D.  TWO  TERM  APPROXIMATION 


In  order  to  proceed  with  the  analysis  so  as  to  derive  a 
signal  processing  algorithm  from  Equation  (2.19),  we  will 
assume  that 

oo  - 

if  I  T  «  1  (2-20> 

o  k=-°° 

In  Appendix  B,  we  will  show  under  what  conditions  this 
assumption  is  valid.  It  turns  out  that  the  assumption  of 
Equation  (2.20)  is  essentially  equivalent  to  an  assumption  of 
low  SNR,  i.e.,  E/Nq  <<  1. 

Since 

ex  -  1  +  x  if  x  <<  1  (2.21) 

provided  Equation  (2.20)  is  valid,  it  is  possible  to  approxi¬ 
mate  the  exponential  term  of  Equation  (2.19)  by  two  terms  only 
so  that  the  test  can  be  written  as 

0  .  H. 

,2vk.  1 

0°  oo  —  J  - \ 

J  (1  +  jp  l  VkRke  T  )f  (X)dX  *  n  exp{E/No)  (2.22) 

-oo  o  k=-«>  „ 

H0 

or  equivalently 


l  W 

k=-oo 


.  2irk. 

*  “3— A 

/  e  f A ( X ) dX 


«] 

> 

< 

Hr 


N 


-y[n  exp{E/NQ}-l]  = 


(2.23) 
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provided  that  Tq  =  nT,  where  n  is  an  integer.  This  result 
also  valid  even  if  n  is  not  an  integer,  provided  that  n  >> 
Thus  applying  this  result  to  Eq.  (2.14)  results  in 

T 

/°v2(  t-X)dt  -  T  I  VkV*  -  T  I  |Vk|2  $  E 

0  k=-°°  k=-°° 


Observe  that  E  is  the  energy  of  the  signal  v(t)  for  0  £  t 
Furthermore , 


T 

o 

/  r(t)v(t-A)dt 

0 


Vke 


r  O  j^-kt 

/  r(t)e  dt 

0 


k=-<» 


VkRke 


-3 


where 


*k 


T  .  2  irk. 

,  °  1 1 

I  r(t)e  dt 

0 


So  the  test  of  Equation  (2.6)  becomes 


00 


/ 


exp{ 


N 


l 

k=-°° 


_  ,  2TTk  ' 

W  T  >£A(X)aX 


*  n  exp{E/NQ} 

% 


is 

1. 


2.16) 

T  . 
o 


2.17) 


2.18) 


2.19) 
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L1(r) 

LftTrl* 


/  exp{— -  /  r{t)v(t-X)dt} 

5  0 


T  "1 


.  exp{-  ~  /  v2(t-X)dt)fA(A)dX  *  n 

o  0  u 


(2.11) 


Since  v(t)  is  T-periodic,  it  can  be  expressed  in  terms  of  an 
exponential  Fourier  series,  namely 


v(t)  =  l  V.  e 


k=-°° 


(2.12) 


where 


T/2  -j^-kt 

i-  /  v  ( t)  e  dt  for  all  k 

T  -T/2 


(2.13) 


Using  the  Fourier  series  expansion,  we  have 


00  00 


/  v  (t-X)dt  = 


I  ,  l  W 

fc=  — 00  £,=-00 


-jTT(k+A)X  To  j^(k+£)t 


(2.14) 


However  note  that 


T  •  2  TT  .  . 

o  D-TB-(k+4)t 

/  e  dt 


if  Z  = 


(2.15) 


0 


otherwise 


The  likelihood  ratio  test  is  therefore 


L1(r) 


H: 

> 

< 

H, 


(2.6) 


where 


and 


P(H0)  C10 

n  '  ^1?"  coi 


P{H.}  i  =  1,0 


(2.7) 


(2.8) 


are  the  prior  probabilities  of  occurrence  of  hypotheses  H^, 
and 


i,j  =  0,1 


(2.9) 


are  the  costs  associated  with  making  decisions  about  hypotheses 
•  For  a  minimum  probability  of  error  receiver 


1  i  ¥  j 

0  i  =  j 


(2.10) 


Taking  advantage  of  some  simplifications  in  writing  the  LR,  it 
is  simple  to  show  that  the  LR  becomes 


likelihood  ratio  (.LRl,  with  a  threshold  n  and  decide  in  favor  of 
hypothesis  if  LR  is  greater  than  n ,  or  in  favor  of  Hg  if 
LR  is  less  than  n. 


Due  to  the  random  variable  X,  v(t-X)  itself  is  a  random 
process.  If  we  condition  on  X,  i.e.,  we  assume  some  fixed 
value  of  X, then  v(t-X)  is  completely  known.  However,  the 
observed  signal  r(t),  conditioned  on  X,  because  of  the  noise 
n(t),  is  a  Gaussian  random  process.  In  other  words  condition¬ 
ing  on  X,  under  either  hypothesis,  r(t)  is  Gaussian  so  that 
the  likelihood  functions  can  be  expressed  in  closed  form  and 
are  given  by  [Ref.  4] 


To 

L  (r)  =  F  exp{-  /  r2(t)dt} 

0  o  0 


(2.3) 


and 


1^  (r | X) 


F  exp 


{- 


T 

O  ? 

/  [r ( t ) -v ( t-X ) ]  at } 

0 


(2.4) 


The  constant  F  is  of  no  consequence  here  as  it  will 
cancel  once  the  LR,  namely  L1(r)/LQ(r)  is  evaluated.  Observe 
that 


oo 

L:(r)  =  /  L1(r|X)fA(X)dX 


(2.5) 


C.  COMPUTATION  OF  THE  LIKELIHOOD  RATIO 


In  this  section,  principles  of  statistical  communication 
theory  are  applied  to  the  derivation  of  a  decision  rule  which 
will  lead  to  the  design  of  a  receiver  that  optimally  detects  (with 
minimum  probability  of  error,  P  )  a  periodic  signal  with  random 
time  delay,  in  the  presence  of  additive  white  Gaussian  noise  (WGN) 
under  the  assumption  of  low  signal  to  noise  ratio  conditions. 

We  start  the  problem  by  considering  two  hypotheses  H-^  and 
Hq  such  that  during  the  observation  interval  (0,Tq),  under 
hypothesis  H^ ,  we  assume  that  the  periodic  signal  with  unknown 
time  delay  is  present  and  under  the  second  hypothesis,  i.e., 

Hq,  there  is  no  signal  present.  In  both  cases,  the  effect  of 
the  WGN  must  be  considered.  Thus  symbolically, 


H^:  r(t)  =  v(t-A)  +  n(t) 

Hq  :  r  ( t )  =  n  ( t) 


0  <  t  <  T 
—  —  o 


(2.2) 


where 


r(t)  =  the  signal  at  the  front  end  of  the  receiver; 

v(t)  =  the  T-periodic  deterministic  signal  with 
A  a  random  variable  modeling  the  unknown 
time  delay  of  the  waveform;  and 

n(t)  =  a  sample  function  of  a  white  Gaussian  process 
with  zero  mean  and  two-sided  power  spectral 
density  level  Nq/2  watts/Hz. 

In  order  to  satisfy  the  minimum  probability  of  error  cri¬ 
terion,  the  optimum  decision  rule  is  obtained  by  comparing  the 
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B.  PROBABILITY  DENSITY  FUNCTION  OF  A 

The  type  of  random  signal  parameter  dealt  with  in  this 
thesis  is  the  phase  or  time  delay  of  an  arbitrary  periodic 
waveform  of  period  T.  Generally,  the  period  T  of  the  signal 
is  very  much  shorter  than  the  time  duration  Tq  of  the  signal 
and  therefore  it  is  extremely  difficult  to  predict  the  time 
delay  at  the  receiver  of  such  a  signal.  In  such  cases  it  is 
reasonable  to  model  the  time  delay  'A',  as  random  variable 
having  some  density  function  that  reflects  the  degree  of 
uncertainty  that  exists  about  A. 

A  useful  p.d.f.  for  the  time  delay  [Ref.  3]  is  given  by 

2tt  A 
mcos— 

£A(A!"°  -  eTl0(m)  X  i  lT/2l  (2'1! 

0  ^  m  £  00 

where  the  independent  parameter  'm'  determines  the  degree  of 
uncertainty  about  the  time  delay  A,  and  IQ(m)  is  a  modified 
Bessel  function  of  the  first  kind  of  order  zero. 

For  m  =  0,  we  see  that  f  ^  ( A }  =  ^  for  A  <_  |  T/2  |  ,  so  that 
A  is  uniformly  distributed  over  the  period  T.  As  m  increases, 
the  uncertainty  of  the  random  variable  A  decreases  and  the 
p.d.f.  approaches  the  shape  of  a  normal  density.  As  m  -*■ 
f^(A)  tends  to  an  impulse  so  that  the  time  delay  uncertainty 
tends  to  zero  and  the  modeled  random  waveform  approaches  that 
of  a  deterministic  or  completely  known  signal.  The  shape  of 
the  p.d.f.  as  a  function  of  m  is  shown  in  Fig.  2.1. 

The  derivation  of  the  mean  and  variance  of  this  density  func 


tion  is  presented  in  Appendix  A. 


II.  NONCOHERENT  RECEIVER  ANALYSIS 

A.  BASIC  CONCEPT 

Detection  of  arbitrary  periodic  signals  with  random  time 
delay  is  a  problem  of  significant  importance  in  radar  and 
communications,  which  has  not  received  a  great  deal  of  attention 
in  the  past.  The  most  relevant  documentation  related  to  this 
problem  appears  in  the  radar  literature,  where  the  problem  is 
formulated  as  noncoherent  detection  of  a  radio  frequency  sine 
wave  burst  in  the  presence  of  noise  [Ref.  1J . 

In  communication  systems,  signals  are  transmitted  to  a 
receiver  via  some  medium  connecting  the  transmitter  to  the 
receiver  called  the  channel.  Transmitted  signals  undergo 
distortion  in  the  channel  as  well  as  in  the  receiver  itself.  In 
many  cases  this  distortion  is  caused  by  physical  processes 
which,  because  of  their  complexity  must  be  modeled  via  statis¬ 
tical  means,  i.e.,  random  variables  and/or  processes.  As 
previously  pointed  out  in  the  receiver  itself,  noise  is 
unavoidably  added  to  the  signals  to  cause  further  distortion 
and  uncertainty.  It  is  often  found  in  practice,  however, 
that  the  uncertainty  created  by  the  noise  in  a  receiver  results 
in  the  introduction  of  signal  uncertainties,  modeled  as  random 
processes  having  known  waveshapes  but  random  parameters  such 
as  random  amplitudes,  frequencies  and/or  phases.  The  problem 
of  detecting  signals  of  random  time  delay  has  been  partially 
considered  in  Ref.  2. 


have  been  carried  out,  and  wherever  possible,  comparisons  are 
made  with  optimum  systems  operating  in  similar  environments. 

The  results  of  Chapters  II  and  III  involving  receiver  Pg 
are  analyzed  and  interpreted  via  the  use  of  tables  and  graphs 
in  Chapter  IV. 

The  fifth  and  final  chapter  presents  some  general  conclu¬ 
sions  to  be  derived  from  the  work  carried  out  in  preparation 
of  this  thesis,  and  some  suggestions  for  future  analysis  in 
this  general  topic  area  are  given. 


detect  an  arbitrarily  shaped  periodic  waveform  with  random 
time  delay,  is  carried  out.  The  random  time  delay  is  assumed 
to  have  some  known  probability  density  function  (p.d.f.). 

Two  types  of  p.d.f . 's  are  assumed  for  the  random  time  delay  A. 
Namely,  a  uniform  p.d.f.  and  a  (.parameter  varying)  non-uniform 
p.d.f.  have  been  considered,  and  their  effect  on  the  receiver 
probability  of  error  (Pg)  has  been  studied. 

The  design  of  optimum  (in  minimum  Pg  sense)  receivers  is 
based  on  the  likelihood  ratio  and  the  assumption  of  low  signal- 
to-noise  ratio  (SNR)  operations.  The  design  of  suboptimum 
receivers  is  based  on  heuristic  approaches  that  intuitively 
yield  reasonably  good  performance. 

Chapter  II  of  this  thesis  presents  statistical  communication 
theoretic  principles  that  have  been  applied  to  the  analysis, 
design,  and  performance  evaluation  of  noncoherent  receivers. 
Approximations  are  made  that  lead  to  a  decision  rule  or  equiva¬ 
lently  a  receiver  structure.  A  justification  and  meaning  of 
such  approximations  is  presented  in  Appendix  B.  It  must  be 
pointed  out,  that  at  high  SNR,  most  receivers  designed  strictly 
from  heuristic  considerations  perform  adequately.  At  low  SNR 
however,  the  receiver  design  must  be  optimum  so  as  to  not 
further  degrade  marginal  operating  conditions.  Chapter  II 
addresses  this  issue. 

Chapter  III  is  devoted  to  the  analysis  and  design  of  sub¬ 
optimum  receivers  where  simple  (heuristic)  techniques  have 
been  applied  to  the  design  process.  Performance  evaluations 


Var{  £  |  H1 ,  A  } 


E{  [  U|HlfX)  -  E{  £  |  H1 ,  X  }  ]  } 


=  E{ [Vn  /  [v ( t-X )  +  n (t) ] dt  -  VQ  /  v (t-X) dt]  } 


o  °  2 

-  v£  /  /  E{n  (t)  n  (t)  }dt  dx  =  -j  VQ  Tq 

=  Var{£|HQ}  =  (2.34) 

Equation  (2.33)  shows  that  the  conditional  mean  m1  is  indepen¬ 
dent  of  X  and  Equation  (2.34)  also  shows  that  the  conditional 
variance  of  £  is  independent  of  X . 

Hence 


fL(£|H1,X)  =  fL(«-|H1) 


X2 

- ( £-m. ) 

-  exp{  - j—} 

/27c£  2o£ 


Also,  from  Equations  (2.32)  and  (2.34) 


fL(*'V 


-  exp{- 


(2.35) 


(2.36! 


The  p.d.f.'s  of  Equations  (2.35)  and  (2.36)  are  plotted  in 
Fig.  2.3. 


Fig.  2.3  Probability  Density  Function  of  i 


From  Equations  (2.35)  and  (2.36)  and  Fig.  2.3,  it  is  clear 
that  the  two  p.d.f.'s  are  completely  identical  with  different 
means.  In  order  to  have  minimum  probability  of  error,  the 
threshold  rf  should  be  the  point  where  the  two  p.d.f.'s  inter¬ 
sect.  That  is. 


n  = 


mi  +mo 


„2m 
Vn  T 
0  o 


(2.37) 


This  threshold  differs  from  that  given  by  Equation  (2.24)  . 
We  must  however  keep  in  mind  that  the  results  were  obtained  by 
approximations  on  the  LR  test  of  Equation  (2.19). 

Receiver  probability  of  error  Pq,  is  thus  given  b.y 


Pe  -  P(H0>Pp  +  P(H1>Pm 


(2.38) 


where 


Probability  of  false  alarm; 


PM  =  Probability  of  miss;  and 


P(H^},  i  =  0,1  are  the  prior  probabilities. 


Since  we  have  assumed  equal  prior  probabilities,  Pg  becomes 
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From  Equations  (2.24),  (2.33)  and  (2.34),  we  get 


=  i{erfc*[ 
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Observe  that  if  v(t)  has  a  very  strong  d.c.  component  in 

2 

comparison  to  its  harmonics,  then  Vqt0/n0  is  approximately 
equal  to  the  SNR  defined  as  E/Nq .  Then  Equation  (2.40)  becomes 
a  function  of  SNR  only.  Otherwise,  it  is  a  function  of  both 
the  SNR  and  the  V^T  /N  ratio. 

ii  n 


2 .  Nonuniform  p.d.f. 

For  this  case,  the  p.d.f.  of  X  as  given  by  Equation 

(2.1)  is 
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It  is  known  that  [Ref.  5] 
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Therefore,  Equation  (2.1)  becomes 
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From  Equation  (2.23)  we  have  the  test 
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or  equivalently 


where 
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Substituting  Equation  (2.42)  into  (2.44),  we  get 
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Note  that  both  integrals  in  Equation  (2.45)  are  zero 


(2.43] 


(2.44) 


(2.45) 


unless 


p  =  ±k.  For  p  =  k,  first  integral  in  brackets  is  T,  and  the 
second  integral  becomes  zero.  For  p  =  -k,  the  first  integral 
is  zero  but  the  second  integral  is  T. 


Therefore 
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(2.46) 


since 


Hence  the  test  of  Equation  (2.43)  becomes 
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(2.47) 


The  test  can  be  simplified  further  by  letting 
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(2.48) 


Notice  that  this  receiver  can  operate  on  signals  both 
with  and  without  d.c.  component.  When  a  signal  has  zero  d.c. 
term,  the  receiver  of  Fig.  2.4  can  be  simplified  by  eliminating 
the  uppermost  branch. 

The  performance  of  this  receiver  is  obtained  by  using 
the  fact  that  conditioned  on  X,  under  either  hypothesis,  g 
is  a  Gaussian  random  variable  with 

E{g  |  Hq}  =  0  (2.51) 

and 
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The  integral  of  Equation  (2.52)  has  been  analyzed  in  Equation 
(2.15).  Therefore 

T  j2TTkx 

■tslHi.U  -  T-THT  !  !vkl2lk<™»e  '  =  hU>  (2‘ 

0  K=-°° 

Furthermore , 


33 


X 

00  Q 

Var{g|H1,A}  =  E{  [  V^On)  /  [v ( t-A )  +n (t)  ]  e 


.  2irk, 
D-rp-t 


dt 


.  27Tk, 


'  W  kl,W">  0I  T  tdt>2> 


oo  To  j  2jrk_ t 

E1(JFsr  k.LvkIk(n)  0J  n(t,e  T  dt)2) 


00  oo  O 

—  l  l  V.v  I  (m)I.(m)  //  Efn(t)n(T)} 

I q  (m)  k=-»  L=-c°  K  *■  K  *  0 


.  2Trk,  -  2tt £ 

D-sr-t  3-st-t 
e  1  e  1  dt  dx 


N  oo  oo  To  jis-(k+ll)t 

4  l  l  Wk(m)Vm)  /  e  dt 


— 2  2 

In  (m)  k=-<»  i=-oo 


T  N 
o  o 

0 T2 ,  ,  .  4  ''k1  ‘k 

2 IQ  (m)  k=-°° 


I 
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With  p  =  P { Hq }  and  1-p  =  PfH^},  we  have  an  expression  for  the 


conditional  P  given  by 
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where  n*/  h(A)  and  a  are  given  by  Equations  (2.24),  (2.53) 

g 

and  (2.54),  respectively. 

The  actual  receiver  P&  is  obtained  by  integrating 

P  (A)  over  the  p.d.f.  of  A,  resulting  in 
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Observe  that 
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Also 
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E.  THREE  TERM  APPROXIMATION 

Whenever  Equation  (2.20)  is  not  satisfied  to  the  degree 
that  Equation  (2.21)  is  valid,  it  is  possible  to  proceed  with 
the  analysis  beyond  Equation  (2.19),  if  in  Equation  (2.19)  we 
expand  the  exponential  term  and  approximate  it  with  three 
terms.  Since 

2 

exp{x}  -  1+x+jp,  if  x  <<1  (2.59) 


in  order  to  have  the  approximation  error  not  exceed  10%,  we 
must  have  x  £  0.79.  Using  the  three  term  approximation  of 
Equation  (2.59)  and  applying  it  to  Equation  (2.19),  the  test 
then  becomes 
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or  equivalently 
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where 
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Here  again,  two  cases  must  be  considered: 

1.  Uniform  p.d.f.  on  X 

2.  Nonuniform  p.d.f.  on  X 
1 .  Uniform  p.d.f. 

For  this  case  the  first  term  on  the  left  hand  side  (L.H.S.) 
of  the  inequality  in  Equation  (2.61)  has  already  been  analyzed 
in  Equations  (2.26)  through  (2.28).  This  first  term  is  just 
VqRq.  The  second  term  of  Equation  (2.61),  i.e., 
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becomes  (with  the  aid  of  Equation  (2.15)) 
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The  test  therefore  becomes 
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or  equivalently  (using  Equations  (2.13)  and  (2.18)) 
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This  test  leads  to  the  receiver  structure  shown  in  Fig.  2.5. 

Note  that  this  receiver  utilizes  both  the  d.c.  and  non- 
d.c.  components  of  the  periodic  signal.  However,  in  practice 
there  are  many  cases  where  the  d.c.  component  of  the  signal  is  zero. 
For  those  cases,  the  test  reduces  to 


4  i  ivki2iv2  : 
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(2.66) 


and  the  corresponding  receiver  can  be  implemented  by  simply 
eliminating  the  upper  most  branch  of  the  jeceiver  in  Fig.  2.5. 
2.  Nonuniform  p.d.f. 


Our  starting  point  here  is  Equation  (2.11),  repeated 


here  for  convenience, 
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where 
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are  the  signal  discrete  cross  correlation  terms.  Thus, 
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Now  conditioned  on  both  and  A,  £  is  also  a  Gaussian 
r.v.  with 
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If  we  consider  the  two  hypotheses,  as  done  in  Chapter  II, 
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then  conditioned  on  hypothesis  Hq  ,  £  is  a  Gaussian  random 
variable  (r.v.)  with 
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Fig.  3.1.  N-Correlator  Receiver 


these  conditions  than  when  the  correlators  are  fed  by  low 
level  noise  only.  Thus,  if  Z  is  greater  than  a  prescribed 
threshold,  then  the  presence  of  a  signal  with  a  random  time 
delay  is  declared;  otherwise  only  the  presence  of  noise  is 
declared.  Here  we  assume  for  convenience  that  the  observation 
time  Tq  is  an  integral  multiple  of  T,  the  period  of  the 
signal  v(t).  Observe  furthermore  that  we  are  not  required  to 
assume  low  SNR  conditions  as  done  in  Chapter  II. 

In  order  to  analyze  this  receiver  in  more  detail,  observe 

that 
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III.  ALTERNATIVE  SUBOPTIMUM  RECEIVERS 


In  the  previous  chapter,  optimum  detectors  based  on 
likelihood  ratio  tests  were  derived  under  conditions  that  are 
equivalent  to  low  SNR  conditions  for  signals  having  random 
time  delay  in  the  presence  of  additive  white  Gaussian  noise. 
The  tests  derived  gave  rise  to  different  receiver  structures 
which,  whenever  possible,  were  analyzed  in  terms  of  their 
error  probability  performance.  Both  the  receiver  structures 
and  their  performance  were  seen  to  be  a  function  of  the  uncer¬ 
tainty  in  the  signal  time  delay  or  more  precisely,  a  function 
of  the  p.d.f.  of  the  signal  time  delay  A. 

In  this  chapter,  two  suboptimum  receivers  are  analyzed. 

The  approach  used  here  is  different  in  the  sense  that  re¬ 
ceivers  are  derived  via  strictly  heuristic  means  and  their 
performances  are  then  analyzed.  These  receivers  are: 

A.  N-correlator  Receiver 

B.  Estimator-correlator  Receiver 

A.  N-CORRELATOR  RECEIVER 

This  receiver  is  shown  in  Fig.  3.1.  The  idea  is  to  corre¬ 
late  the  received  signal  r(t)  with  N  delayed  reference  signals 
generated  at  the  receiver.  If  a  signal  is  present,  presumably 
at  least  one  correlator  branch  would  produce  a  large  output  so 
that  the  sum  l  could  be  guaranteed  to  be  much  larger  under 
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Observe  that  Equation  (2.90)  (as  well  as  Equations 
(2.57)  and  (2.58))  involve  summations  whose  index  runs  from 
1  to  infinity.  In  practice,  it  is  not  possible  to  compute 
infinite  sums.  However,  due  to  the  fact  that  for  all  reason¬ 
able  signals  v(t) ,  the  magnitude  of  the  coefficients  gets 
smaller  as  k  gets  bigger,  the  infinite  sum  can  be  truncated 
without  introducing  significant  computational  error.  Further¬ 
more,  terms  involving  the  magnitude  squared  of  these  coeffi¬ 
cients  will  tend  to  decay  more  rapidly  allowing  early  truncation 
of  the  infinite  sums. 

G.  DISCUSSION 

In  order  to  carry  out  the  evaluation  of  Equation  (2.90) 

to  a  reasonable  degree  of  accuracy,  the  mathematical  expression 

given  by  Equation  (2.90)  was  evaluated  by  computer  for  fixed 

SNR  values.  Computation  of  Pe  was  carried  out  for  a  fixed 

value  on  the  upper  limit  on  the  index  k.  Recomputation  of 

Pe  was  carried  out  everytime  this  upper  summation  limit  was 

incremented  by  1.  Recomputation  of  Pe  was  stopped  after  an 

incrementation  of  the  upper  limit  in  the  sum  terms  did  not 

yield  an  appreciably  different  value  for  Pg  in  comparison  to 

the  P  value  before  incrementation, 
e 

The  result  of  these  evaluations  has  been  presented  in 
Chapter  IV  in  the  form  of  tables  and  graphs  while  Pg  has  been 
plotted  as  a  function  of  SNR  for  different  values  of  the 
parameter  m.  Recall  that  m  specifies  our  relative  prior 
knowledge  about  the  signal's  random  time  delay. 
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Hence  Equation  (2.57)  becomes 


n* 


E/N 

(ne  °-l) 


4E/N° 

A^I q (m) 


(2.88) 


and  Equation  (2.58)  becomes 


n *-h  (X ) 


N 

(ne  °-l)  - 


4E/N 


(2.89) 


Thus  Equation  (2.56)  for  the  receiver  probability  of  error 
becomes 


Pe  =  p  erfc  * 


E/N 

(ne  °-l) 


4E/No  k£1lvkl2lk<”> 


2  2 
IQ{m)A 


1/2 


+  (l-p)  /  erf* 
0 


f  E 

N  4E/N  °° 

(ne  °-l)-( - 2_y)  (  l  |  V  | 

I n(m) r  k=l 

/  u 

2 _  ,  .  2iTk,  . 

Ik  (m)  cos  T  X) 

[r 

1/2 

2  2 

me  os 


2ttX 


TIq (m) 


dA 


(2.90) 
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j*  .■  :■ 


'!iFj T ! ■ . v  -vj.  w  w.*  "  .■  v  v  v  v 


Note  that 


T/2  „ 

f  v  (t)dt 
-T/2 


a2t 


T  l  lvk! 2 

k=-oo 


(2.84) 


which  implies  that 


I 

k=-co 
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k 


2 


A 


2 


(2.85) 


Furthermore 


I  |Vk|2Ik(m)  =  |V0|2I0(m)  +  2  I  |vk|2Ik(m) 

k=-oo  K—  X 


(2.86) 


so  that 


l  lVkl  Ik(m)e 
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-i^x 


vol  1o<*’ 


+  2  T  |Vk|2Ik(m)  cos^X 


(2.87) 


since 

|vk  1  -  lv-kl;  Jr'1"’  -  W0 
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Fig.  2.7  Periodic  Square  Wave 


The  Fourier  coefficients  of  v(t)  are 


,  T/2 

i[  /  Ae  dt 

1  0 


T/2 


_  •  27rk. 

T  -3-sr-t 

/  Ae  dt] 


irk  . 


Observe  that  V  and  all  the  even  coefficients  are  zero,  i 


V_,  =  0  for  all  k 


and 


V 


±k 


(— )  2 


t 


k  odd 


(2.81) 


e. , 

(2.82) 


(2.83) 
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o  O'* 


(2.79 


n  *r  h  (,X  l  _ 
a 


1/2 

I7TO-<2e/no> 


Thus  the  expression  for  Pe  becomes 


p  erfc* 


(ne 

I1(m) 

T  (mY 

0 


(2E/Nq) 


(1-p)  T/2 

+  TI  (m)  /  erf* 

u  -T/2 
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E/N 


Ix  (m) 
I  o  (m) 


9 T7  1 1  (n\) 


(2E/Nq) 


2tt 

mcos-^ 

(2.80 


Computation  of  Equation  (2.80)  can  be  carried  out  as 
a  function  of  the  SNR  E/Nq,  p  and  m.  This  has  been  carried  out 
and  the  results  are  presented  graphically  in  Chapter  IV.  For 
the  sake  of  computational  simplicity, equal  prior  probabilities 
were  assumed,  i.e.. 


P  = 


1-p  = 


1 

2 


resulting  in  n  =  1 . 

2 .  Square  Wave 

For  the  second  case,  let  v(t)  be  a  periodic  square 
wave  as  shown  in  Fig.  2.7. 


1.  Sine  Wave 


2.  Square  Wave 
1.  Sine  Wave 

For  the  first  case  let 


v  ( t ) 


cos  ^-t 


.  2tt  . 

A  ^  Tt 

J  e  + 


(2.75) 


This  periodic  signal  has  only  two  discrete  components. 
In  other  words,  it  can  be  represented  in  terms  of  its  two 
Fourier  coefficients,  i.e., 


V, 


k  =  ±1 

otherwise 


Thus , 


(2.76) 


I 

k=-°° 


(|)2  + 


(|)2 


(2.77) 


and  from  Equations  (2.57)  and  (2.58),  we  obtain 


n* 

o~ 

g 


-i) 


I-i  (m) 

T^-nzrt2 e/n  ) 


1/2 


(2.78) 


which  is  a  function  of  (t-x)  only.  This  implies  that  for  the 
uniform  p.d.f.  on  A,  the  impulse  response  of  the  filter  turns 
out  to  be  time  invariant.  That  is,  it  depends  only  on  the 
time  difference  (t-t) .  Furthermore,  mY(t)  being  a  constant 
can  be  eliminated  from  the  upper  branch  of  the  receiver  shown 
in  Fig.  2.6. 

F.  ANALYSIS  OF  TWO  SPECIFIC  SIGNAL  WAVESHAPES 

It  was  pointed  out  earlier  that  the  receiver  derived  using 
a  two  term  approximation  on  the  exponential  appearing  in  Equa¬ 
tion  (2.19),  is  unable  to  discriminate  between  the  signal 
present  versus  the  signal  absent  case  for  signals  having  zero 
d.c.  component  when  the  random  time  delay  obeys  a  uniform 
distribution.  When  the  random  time  delay  obeys  a  nonuniform 
p.d.f.,  a  receiver  was  derived  (see  Fig.  2.4)  that  could 
detect  the  presence  of  signals  having  zero  d.c.  component. 

In  this  section,  assuming  a  nonuniform  p.d.f.  on  A,  the 
receiver  performance  will  be  analyzed  for  two  different  signals 
which  are  found  quite  often  in  practice  and  that  have  a  zero 
d.c.  component.  The  two  signals  are 


Fig.  2.6  Noncoherent  Receiver  (Nonuniform  p.d.f .  , 
3  Term  Approximation) 


T/2 

nu,  ( t)  =  ^  /  v(t-A)dA 

T  t-T/2 


jllJit 

e  T  / 

-T/2 


T/2 


dA 


(2.73) 


which  is  a  constant,  and  the  autocorrelation  function  of  y(t) 
given  by  Equation  (2.71)  turns  out  to  be 


Note  that  v(t)  is  a  deterministic  signal  but  as  mentioned 
earlier,  due  to  the  r.v.  A, 

A 

y(t)  =  v(t-A)  (2.69) 

itself  is  a  random  process  with  expected  value 


00 

E{y  (t)  }  =  /  v(t-A)fAU)dA  4  m  (t) 

—  CO 


(2.70) 


and  autocorrelation  function 


E{y (t) y (x )  }  =  /  v(t-A)v(x-A)fA(A)dA  =  Ry(t,x)  (2.71) 

Substitution  of  Equations  (2.70)  and  (2.71) 
in  Equation  (2.68)  leads  to  the  following  test  that  a  receiver 
has  to  perform 

T  T  T  H1 

2  °  n  ^  ^  > 

5-  /  r(t)mv(t)dt  +  -4  /  r (t) [  /  r (x)Rv(t,x) ]dt  „  y  (2.72) 

o  0  Y  0  0  Y  * 

°  H0 

and  the  implementation  of  this  receiver  is  shown  in  Fig.  2.6. 
Here  the  autocorrelation  function  R^(t,x)  can  be  considered 
as  the  impulse  response  of  a  filter  that  has  to  be  designed 
for  detection  of  the  signal. 

Furthermore,  for  the  uniform  p.d.f.  on  A,  the  expected 


value  of  y(t)  given  by  Equation  (2.70)  becomes 


n 


/  exp{tr-  /  r(t)v(t-X)dt}exp{-  /  v  (t-A) dt}f . (A ) dA  * 

5  wo  0  o  0 

H , 


where  n  is  given  by  Equation  (2.7).  Note  that  the  second 
exponential  of  the  inequality  given  by  Equation  (2.11)  has 
already  been  analyzed  in  Equations  (2.14)  through  (2.16). 
Furthermore,  if  we  expand  the  first  exponential  in  the  above 
inequality  and  approximate  it  with  three  terms  as  given  by 
Equation  (2.59),  then  the  test  of  Equation  (2.11)  becomes 


°°  ~  o 

J  J  r  (t)  v(t-A)dt  f.  (A)dA 

-«  o  0  A 


OO  Q 

+  |  /  l~  I  r(t)v(t-A)dt]2fA(A)dA 

-00  o  0 


(2.67) 


where  y  is  given  by  Equation  (2.62). 

Interchanging  the  order  of  integration  in  Equation 
(2.67),  we  obtain 


/  r  (t )  /  v(t-A)f  (A)dAdt 


+  — 2  If  r  ( t )  r  ( x )  /  v(t-A)  v(t-A)  f  ( A )  dA  dtdx  ^  y 


> 
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W 


(2.68) 


(3.10) 


(3.14) 


fL  ( i  |  H1 ,  A )  = 


-u  -  I  an(A))2/2 a2 


l_  e  n-i 


yfi- iroT 


with 


fL  ( £  |  =  /  f  ( «.  1  Hx  ,  A )  f  A  ( X )  dA 


(3.15) 


The  detector  of  Fig.  3.1  performs  the  test 


A  <  Y 


(3.16) 


so  that  its  error  probability,  Pg  is  given  by 


Pe  =  P  (HQ }  /  p  (4|H0)d£  +  P{H.}  /  p  U|H ,)  d£ 

Y  —  oo 


(3.17) 


Let  p  =  P{HQ}  and  1-p  =  P{H^}  so  that  using  equations 
(3.13)  and  (3.15),  equation  (3.17)  becomes 


-U-  l  a  (A)) 

n=l 


00  -i  -l2 /2o2  y  oo  . 

p  /  -z=-e  Wa-p)  /[  /  — z 


y  /2iTa( 


/27  a, 


f A ( A ) dA ] dfl,  (3.18) 
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Interchanging  the  order  of  integration  in  the  second 


(double)  integral  of  equation  (3.18),  we  obtain 
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/  fA<*> 


-U-  l  o  IX)) 
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1 - 
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2c 


dl 


-00  /2ttc, 


dA 


•  N 


/  ffl(A) 
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dA 
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l  an  ( A  ) 

/  f .  ( A )  erf  *|  — -  I  dA 


Thus,  the  receiver  probability  of  error  becomes 


N 


p  =  p  erfc*  +  ( 1— p )  /  f.(A)erfi 


fy-  I.anU) 


n=l 


(3.19) 


dA  (3.20) 


The  threshold  y  can  be  optimized  to  minimize  Pg.  Such  an 

approach  unfortunately  leads  to  a  complicated  integral  equation 

that  is  not  easily  solvable,  in  practice  however,  P  can  be 

e 

calculated  on  the  computer  for  y  =  0,  and  then  incrementing 
| y |  until  an  optimum  value  of  y  is  found  for  a  given  SNR. 

It  must  be  pointed  out  that  the  receiver  of  Fig.  3.1  will 
not  always  yield  desirable  results.  For  example,  if  v(t)  is 


sinusoidal,  in  the  absence  of  noise,  l  =  0  which  means 
that  with  noise  present,  the  receiver  is  unable  to  detect 
the  presence  of  the  signal.  Fortunately  for  v(t)  sinusoidal, 
a  better  receiver  is  available  and  its  performance  is  well 
known  [Ref.  4]. 

1 .  Special  Case  of  Triangular  Wave 

Consider  now  a  periodic  waveform  of  triangular  pulses 
and  its  two  arbitrarily  delayed  versions  as  shown  in  Fig. 

3.2.  If  we  represent  one  triangular  pulse  by  p(t),then  we 
can  express  v(t)  as 

M 

v ( t )  =  l  p(t-jT)  (3.21) 

j=0 

In  order  to  find  the  cross  correlation  terms  A  m 

nm 

required  for  the  computation  of  Pg  [equation  (3.19)], 

instead  of  taking  all  the  pulses  in  equation  (3.21),  we 

simply  consider  the  cross  correlation  of  a  single  pulse  of 

period  T  with  a  delayed  version  of  itself  and  then  multiply 

the  result  by  the  total  number  of  pulses  M  within  the  interval 

T  ,  since  v(t)  is  assumed  to  be  periodic  with  MT  =  T  . 
o  o 

In  order  to  make  computation  simpler.  Fig.  3.3  shows  p(t), 
p(t-a)  and  p(t-B),  where  a  <  6  has  been  assumed. 


Fig.  3 


3  One  Period  Restricted  Triangular  Wave  and 
Two  of  Its  Delayed  Versions 
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Observe  that  for  a  <  3 


/  v(t-a) v(t-8)dt  =  M 
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equation  (3.22)  becomes 
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n  nm  , 

2N  '  ^ 


m 


2N  +  3) 


A 


nm 


(3.23) 


provided  n  <  m.  For  n  =  m  or  equivalently  a  =  8,  the  quantity 
inside  the  parenthesis  in  equation  (3.23)  becomes  j. 

Similarly  for  a  >  8,  equation  (3.22)  becomes 

T 

/  v  ( t-a)  v  ( t-8 )  dt  =  ma2(^^-^ - a6T  +~-  +^-)  (3.24) 
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provided  n  >  m.  Here  again  for  n  =  m,  the  quantity  inside 
the  parenthesis  in  equation  (3.24)  becomes  y. 

To  proceed  further  with  the  analysis,  the  double 
summation  of  equation  (3.8)  is  broken  into  three  summations 
to  account  for  n  <  m,  n  =  m  and  n  >  m.  Thus 


N 


l 

n=l 


N 


l 

m=l 


N  N-l  N  N  N-l 

l  l  Anm  +  l  Ann  +  2  £ 

m=n+l  n=l  nm  n=l  nn  n=m+l  m=l 


A 

nm 


(3.26) 


Observe  that  the  L.H.S.  of  equation  (3.26)  represents 
the  sum  of  the  elements  of  an  N  «n  matrix.  Whereas  the  middle 
term  on  the  R.H.S.  of  this  equation  denotes  the  sum  of  the 
diagonal  elements  of  the  N  xN  matrix,  the  first  term  is  the 
sum  of  all  the  elements  in  the  upper  triangular  part  and  the 
third  term  is  the  sum  of  all  the  elements  in  the  lower  triangu¬ 
lar  part  of  the  N  *N  matrix. 

In  order  to  evaluate  the  first  term  on  the  R.H.S.  of 
equation  (3.26),  use  of  equation  (3.23)  (n  <  m  case)  yields 
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(3.27) 


The  second  term  on  the  R.H.S.  of  equation  (3.26) 


becomes 
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l,  A„n 
n=l 


,2T3  i  I 


ma2t3  [|] 


(3.28) 


since  A  „  =  ■=■  for  n  =  m. 
mn  3 

Finally,  substitution  of  equation  (3.25)  (n  >  m  case) 
in  the  third  term  on  the  R.H.S.  of  equation  (3.26)  yields 


N  N-l 
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(3.29) 


Adding  the  results  of  equations  (3.27),  (3.28)  and 


(3.29)  yields 


N  N 
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(3.30) 


Thus  from  equations  (3.8)  and  (3.30),  the  conditional 


variance  of  £  becomes 
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since 


Tq  >P  ^2 
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is  the  energy  of  v(t)  for  0  £  t  £  Tq . 

In  order  to  compute  the  probability  of  error  P0 ,  we 
N 

need  to  evaluate  7  a  (A)  given  by  equations  (3.9)  and  (3.11) 

n=l  n 


N  No  nT 
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i=l  n  n=l  0 


(3.33 
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•  .  nT 

A  Fourier  series  expansion  of  v(t-A)  and  v(t  ) 
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yields 
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n=l 


Observe  that  for  k  =  ±iN,  where  i  is  an  integer,  the 
second  summation  involving  the  index  n  yields  N,  otherwise  it 
yields  zero. 

Thus,  equation  (3.33)  can  be  written  as 
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nT  inT 

In  equation  (3.6),  if  we  represent  v(t  — ^-)  and  v(t  — — ) 
in  terms  of  Fourier  series,  then 
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since 
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sin (N0/2)  „  j (N+l ) 9/2 
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(3.37) 


Observe  that  the  R.H.S.  of  equation  (3.36)  is  always  zero 
except  for  k  =  ±iN  where  i  is  an  integer.  Also  observe  that 
for  k  =  ±iN,  the  R.H.S.  of  equation  (3.36)  becomes  an  indeter¬ 
minate  form.  Due  to  that  reason,  differentiating  sin(kfT)  and 

k  7T 

sin(— )  w.r.t.  k  and  then  taking  the  limit  k  -+  iN,  equation 
(3.36)  becomes 
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n-1  m-1  Anm 


T  N2  l  1 V  |  2 
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(3.38) 
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so  that  equation  (3.8)  can  be  written  as 


°l  -  ^TTon2(V0  +  2  lVkN|2) 


(3.39) 


Equations  (3.35)  and  (3.39)  are  of  general  form  and  can 
be  applied  to  any  signal,  not  just  a  triangular  waveform. 

The  Fourier  coefficients  of  the  periodic  triangular 

wave  are 


^  k  =  0 


otherwise 


(3.40) 
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4tt  k 


k  =  0 


otherwise 


(3.4i: 


so  that  the  use  of  equations  (3.31)  and  (3.35)  yields 
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(3.42) 


where  y*  =  y/N 
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1  +  1/3N* 


(3.45) 


Finally,  the  expression  for  P  ,  after  substituting  equations 
(2.1),  (3.42)  and  (3.44)  in  equation  (3.20),  becomes 
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Figure  4.1  shows  the  performance  of  this  receiver  as 
a  function  of  SNR.  This  graph  was  obtained  through  numerical 
evaluation  of  Equation  (2.40)  under  the  assumption  that  v(t) 
has  a  strong  DC  component  in  comparison  to  its  harmonics, 
equal  prior  probabilities,  and  uniform  p.d.f.  on  the  r.v.  X.  Due 
to  the  constraint  imposed  by  Equation  (.2.20),  Equation  (2.40) 
is  valid  up  to  an  SNR  of  about  -4.1  dB .  Thus  Fig.  4.1  is 
plotted  only  up  to  the  SNR  of  -4.1  dB . . 

Observe  that  Pe  is  at  a  very  high  level  of  about  0.3. 
However,  it  decreases  as  SNR  increases,  as  expected.  While 
it  is  clear  that  no  receiver  could  operate  with  an  error 
probability  of  0.3,  the  figure  does  give  an  indication  of  the 
performance  level  that  can  be  expected  at  low  SNR's.  For 
high  SNR's,  it  is  expected  that  P  will  continue  to  decrease 
to  reasonable  levels.  However,  for  high  SNR's,  the  receiver 
of  Fig.  2.2  is  no  longer  optimum. 

2 .  Receiver  Operating  on  Arbitrary  Signals  (Nonuniform 
p.d.f.  on  X) 

The  two-term  approximation  on  the  exponential  appearing 
in  Equation  (2.19)  and  application  of  the  nonuniform  p.d.f.  on 
X  given  by  Equation  (2.1)  resulted  in  the  receiver  of  Fig.  2.4 
that  could  operate  on  signals  having  either  zero  or  non-zero 
DC  component.  Observe  that  as  m  -*■  0 ,  the  receiver  of  Fig. 

2.4  becomes  equal  to  the  receiver  of  Fig.  2.2.  Two  specific 
signals,  a  sine  wave  and  square  wave,  both  having  zero  DC 
component,  were  analyzed. 
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IV.  DISCUSSION  OF  GRAPHICAL  RESULTS 


This  chapter  presents  graphical  results  obtained  by 
applying  the  analytical  results  of  the  previous  chapters  to 
specific  examples  and  carrying  out  the  required  computations 
on  the  computer.  The  plots  are  intended  to  display  receiver 
performance  in  terms  of  probability  of  error  (Pq)  as  a  function 
of  SNR.  Some  curves  of  Pg  versus  comparator  threshold  level 
have  also  been  included  in  order  to  show  the  dependency  of  Pe 
on  the  threshold  which  in  turn  depends  on  SNR.  The  graphical 
results  displaying  performance  of  the  optimum  receivers  are 
presented  first,  followed  by  graphical  results  displaying 
performance  of  the  suboptimum  receivers  analyzed. 

A.  GRAPHICAL  RESULTS  FOR  OPTIMUM  RECEIVERS 

1 .  Receivers  Operating  on  Signals  with  Nonzero  DC 
Component  (Uniform  p.d.f.  on  A) 

In  Chapter  II,  using  statistical  communication  theoretic 
principles,  a  simple  integrator  receiver  (Fig.  2.2)  was  derived 
that  could  discriminate  between  signal  plus  noise  and  noise 
only  hypotheses  provided  the  signal  had  a  non-zero  DC  component. 
Recall  that  the  receiver  derivation  and  the  subsequent  per¬ 
formance  evaluation  carried  out  was  made  possible  by  the 
two-term  approximation  on  the  exponential  appearing  in  Equation 
(2.19).  Equivalently,  low  SNR  conditions  were  assumed  so  that 

the  P  results  are  only  valid  for  low  SNR's, 
e 
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summation  in  equations  (3.35)  and  (3.61)  turned  out  to  be  of 
such  a  form  that  it  could  be  expressed  in  closed  form.  . 

For  the  sine  wave,  summation  disappeared  due  to  the  fact 
that  its  coefficients  exist  only  for  k  =  ±1.  These  results 
in  the  form  of  graphs  are  presented  and  discussed  in  Chapter 


so  that  equation  (3.59)  becomes 
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Observe  that  equations  (3.35)  and  (3.62)  involve 
summations  whose  indices  run  from  1  to  infinity.  In  practice, 
it  is  not  possible  to  compute  infinite  sums.  However,  due  to 
the  fact  that  for  all  reasonable  signals  v(t),  the  magnitude 
of  the  coefficients  gets  smaller  as  k  gets  bigger,  this 
infinite  sum  can  be  truncated  without  introducing  significant 
computational  error. 

Fortunately  for  triangular  and  square  wave  signals 
considered  as  special  cases  in  this  chapter,  the  argument  of 
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Thus  equation  (3.60)  becomes 
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and  for  p  =  (1-p)  =  j,  equation  (3.58)  becomes 
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As  mentioned  earlier  in  Chapter  II,  m  =  0  corresponds 
to  a  uniform  p.d.f.  on  A,  i.e.  p(A)  =  so  that  the  error 
probability  expression  for  m  =  0  becomes 
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2 .  Square  Wave 

Consider  now  the  second  special  case  in  which  v(t) 
is  a  periodic  square  wave  as  shown  in  Fig.  2.6.  In  Chapter  II 
we  showed  that 


V, 
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(2A)2 

'  kfT 


for  k  =  odd 


for 


k  =  even 


As  before,  determination  of  the  optimum  threshold  by  solution 
dpe 

of  =  0  involves  solving  an  equation  of  a  form  similar  to 
that  of  equation  (3.47). 

In  general,  solutions  in  closed  form  do  not  appear  tracta¬ 
ble.  A  search  for  an  optimum  y  may  be  performed  or  solutions 
for  specific  cases  may  be  possible.  This  is  illustrated  by 
considering  two  special  cases  in  which  v(.t)  is  a 

1.  sine  wave 

2 .  square  wave . 

1 .  Sine  Wave 

As  a  first  special  case,  let 
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so  that 
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Thus  the  expression  for  Pg  is  given  by 
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where 
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as  before, 


After  a  change  of  variables,  equation  (3.57)  can  be 
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It  can  be  shown  without  a  great  deal  of  difficulty  that 
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where  V  are  the  exponential  Fourier  series  coefficients  of 
v(t),  and  E  is  the  energy  of  v(t)  for  0  <  t  <  T. 


Delay 
TqScc . 


Noncoherent  Correlator  Receiver 


Clearly,  the  computer  could  be  used  to  search  for  a  value  of 
y  that  satisfies  the  equality  prescribed  in  equation  (3.48). 

The  approach  that  has  been  used  in  trying  to  obtain 
the  optimum  threshold  involves  using  the  computer  to  evaluate 
equation  (3.46)  as  a  function  of  y  for  fixed  SNR.  Thus  the 
value  of  y  for  which  Pe  is  minimum  yields  the  optimum  thres¬ 
hold  for  that  particular  value  of  SNR.  This  procedure  is 
repeated  for  different  values  of  SNR.  The  resulting  values 
are  then  used  to  plot  Pg  vs.  SNR.  These  curves  have  been 
plotted  using  optimum  threshold  settings  and  the  results  have 
been  summarized  and  discussed  in  Chapter  IV.  (Observe  that 
the  threshold  y  is  SNR  dependent.  This  situation  is  typical 
of  noncoherent  signal  detection  problems.) 

B .  ESTIMATOR-CORRELATOR  RECEIVER 

The  second  suboptimum  receiver  being  proposed  is  the  re¬ 
ceiver  structure  shown  in  Fig.  3.4.  It  is  basically  an 
estimator-correlator  receiver  in  which  first  a  coarse  estimate 
of  the  time  delay  in  v(t)  is  made  and  then  this  information 
is  used  to  process  the  delayed  signal  r(t)  coherently  with 
an  estimated  reference. 

In  this  receiver  the  largest  output  of  the  N  correlator 
branches  is  used  as  an  estimate  of  the  signal  phase  by  the 
logic  choice.  This  estimate  is  used  to  provide  an  accurate 
local  reference  with  which  to  perform  the  correlation  detection 
operation.  The  output  £  of  the  main  correlator  is  used  for 
threshold  comparison  in  order  to  make  the  decisions. 
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As  mentioned  earlier,  the  receiver  of  Fig.  3.1  has 

been  designed  via  ad-hoc  methods,  so  that  the  optimum  threshold 

setting  is  not  readily  obtainable.  The  difficulty  in  applying 

an  analytical  approach  to  the  problem  of  finding  the  optimum 
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threshold  arises  in  the  solution  of  -  0  for  y.  Since 
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Figure  4.2,  obtained  through  numerical  evaluation 
of  Equation  (2.80),  shows  the  performance  of  this  receiver 
as  a  function  of  SNR  when  v(t)  is  a  sine  wave.  Three  differ¬ 
ent  values  of  the  independent  parameter  m  were  considered. 
Recall  that  m  controls  the  level  of  uncertainty  about  the  time 
delay  A.  As  expected,  as  m  increases,  Pe  decreases  for  a 
given  SNR.  For  example,  ai  SNR  of  -5.2  dB  is  required  in 

order  to  achieve  aP  of  0.4  at  m  =  2.  With  m  =  90,  we  need 

e 

an  SNR  of  -9  dB  to  achieve  the  same  Pe .  Note  that  the  m  =  0 
case  has  not  been  considered  here  due  to  the  fact  that  m  =  0 
corresponds  to  a  uniform  p.d.f.  on  A.  The  receiver  for  that 
case  (Fig.  2.4)  is  unable  to  detect  signals  having  zero  DC 
component,  resulting  in  Pe  equal  to  0.5. 

Curves  shown  in  Fig.  4.3  are  the  result  of  numerical 
evaluation  of  Equation  (2.90)  in  which  three  different  values 
of  the  independent  parameter  m  have  been  considered.  The 
m  =  0  case  has  not  been  considered  here  either  due  to  the 
same  reasons  presented  for  the  sine  wave  case. 

Observe  that  Equation  (2.90)  involves  the  computation 
of  infinite  sums  which  can  not  be  carried  out  in  practice. 
However,  due  to  the  fact  that  for  all  finite  power  signals 
v(t),  the  magnitude  of  the  coefficients  gets  smaller  as  k 
gets  bigger  and  thus,  truncation  is  possible.  Inclusion  of 
more  coefficients  in  Equation  (2.90)  corresponds  to  increased 
accuracy  in  the  results.  Furthermore,  Pg  changes  as  k  is 
incremented,  however  up  to  a  point  only  as  further  increase  in 


k  beyond  this  limit  do  not  result  in  an  appreciable  change 
in  PQ.  This  phenomenon  is  illustrated  in  Table  I. 


table  i 

PROBABILITY  OF 

ERROR  VS  »  OF 

COFFICIENTS  SUMMED 

( 

SNR  =  -10  OB) 

*  OF  COFFICIENTS 

probability  of 

SUMMEO 

ERROR 

I 

0.445t  24370+00 

2 

0.444339470+00 

3 

0 .444325590+00 

A 

0.444325430+00 

5 

0.444325420+00 

6 

0. 444325420+00 

7 

0.444325420+00 

a 

0 . 4443 2542 D+ 00 

9 

0 .444  32542  D+  00 

to 

0.444325420+00 

We  see  that  Pq  keeps  on  decreasing  until  the  number  of 
coefficients  summed  is  5.  For  k  >  5,  Pg  remains  unchanged. 
Thus,  for  this  example  the  computation  of  Equation  (2.90)  was 
carried  out  with  sums  truncated  to  5  terms.  The  performance 
curves  for  the  square  wave  signal  are  shown  in  Fig.  4.3 
for  low  SNR  conditions.  With  m  =  2,  an  SNR  of  -4.1  dB  is 
required  in  order  to  achieve  a  Pg  of  0.4  while  for  the  same 
Pe,  with  m  =  90,  the  SNR  required  is  -8.8  dB. 

B.  GRAPHICAL  RESULTS  FOR  SUBOPTIMUM  RECEIVERS 

In  Chapter  III,  two  suboptimum  receivers  were  analyzed. 

The  design  approach  used  in  that  chapter  was  different  in  the 
sense  that  receivers  were  derived  via  strictly  heuristic 


means.  Therefore,  the  optimum  threshold  setting  for  the 
receiver  was  not  readily  obtainable.  However,  a  computational 
approach  was  used  to  solve  the  problem  of  properly  setting 
the  threshold  prior  to  evaluating  Pg  as  a  function  of  SNR. 

1 .  N-Correlator  Receiver 

The  first  suboptimum  receiver  analyzed  in  Chapter  III 
was  the  N-Correlator  receiver  of  Fig.  3.1,  and  its  performance 
with  v(t)  a  periodic  triangular  pulse  was  evaluated.  The 
graphical  results  corresponding  to  the  numerical  evaluation 
of  Equation  (3.46)  are  shown  in  Figs.  4.4  through  4.12.  As 
mentioned  earlier,  using  the  computer.  Equation  (3.46)  has 
been  evaluated  as  a  function  of  the  threshold  y  for  fixed  SNR 
values  starting  at  -10  dB  and  increasing  in  increments  of  5 
dB.  For  a  given  SNR,  the  threshold  value  for  which  Pg  is 
minimum  corresponds  to  the  optimum  threshold  for  that  particu¬ 
lar  SNR.  These  minimum  points  corresponding  to  each  SNR  have 
been  used  as  optimum  threshold  values  in  Equation  (3.46)  and 
the  resulting  P  has  been  plotted  in  Fig.  4.12  as  P  vs  SNR. 

Observe  that  Equation  (3.46)  involves  the  previously 
encountered  parameter  m  and  N,  namely  the  number  of  correlators 
used  by  the  receiver.  Two  different  values  of  each  parameter 
m  =  0,  90  and  N  =  2,  16  have  been  used  in  the  computations 
in  order  to  evaluate  their  effect  on  receiver  performance.  An 
inspection  of  Fig.  4.4  (m  =  0 ,  N  =  16)  and  Fig.  4.6  (m  =  30, 

N  =  16)  reveal  that  for  a  given  N,  m  has  no  significant  effect 
on  P  .  On  the  other  hand  Fig.  4.4  (m  =  0,  N  =  16)  and  Fig. 

4.8  (m  =  0,  N  =  2)  reveal  a  small  effect  of  N  on  Pg  for  a 
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given  value  of  m.  The  actual  effect  of  N  on  P  is  shown  in 

e 

Fig.  4.12  where  Pe  has  been  plotted  as  a  function  of  SNR  for 
two  values  of  N.  Observe  that  the  higher  value  of  N  yields 
better  performance.  To  achieve  a  P£  of  10~^  with  N  =  2,  the 
SNR  required  is  ~  13.5  dB  whereas  for  the  same  Pg  with  N  =  16 
we  need  an  SNR  ~  13.2  dB.  This  result  is  once  again  expected 
as  using  few  correlators  might  result  in  a  small  output  on 
the  correlator  branches  when  v(t)  received  in  r(t)  has  a 
delay  very  much  different  than  the  locally  generated  v(t). 

While  it  is  possible  to  "see"  the  reason  for  the 
independence  on  m  in  Equation  (3.46),  a  heuristic  explanation 
of  this  observed  result  follows  the  following  argument.  The 
receiver  of  Fig.  3.1  correlates  the  incoming  signals  with  local 
replicas  of  v(t)  delayed  by  multiples  of  T/N.  Depending  on 
the  actual  delay  in  v(t),  upon  reception  of  r(t),  one  of  the 
correlator  branches  will  produce  the  largest  output  (or  there 
may  be  a  tie  between  two  branches) .  The  output  produced, 
however,  namely  which  correlator  branch  is  largest,  should 
be  independent  of  the  level  of  uncertainty  about  A.  Clearly, 
as  m  increases,  it  is  easier  to  predict  which  correlator 
branch  output  will  be  largest.  But  this  is  inconsequential 
as  the  receiver  is  only  trying  to  detect  the  presence  of  v(t). 

2 .  Estimator-Correlator  Receiver 

The  second  suboptimum  receiver  analyzed  in  Chapter  III 
was  the  Estimator-Correlator  of  Fig.  3.4.  Its  performance 
for  v(t)  either  a  periodic  sine  wave  or  a  square  wave  was 
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Figure  4.4  Pe  vs  Thld  (m  =  0,  N  =  16)  ,  (SNR  =  “10,  — 
0,5  dB) 
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PE  VSTHLD 

(  1M0.  N»18 ) 


Figure  4.9  Pe  vs  Thld  (m  =  0,  N  =  2) ,  (SNR  =  10,  15, 
16.5,  17  dB) 
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PEVSTHLD 

(  M=30.  N=  2  ) 


10  Pe  vs  Thld  (m  =  30,  N  —  2) ,  (SNR  = 
0,  5  dB) 
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PE  VSTHLD 

( 11=30.  N=2 ) 


Figure  4.11  Pe  vs  Thld  (m 
16.5,  17  dB) 


considered  as  a  function  of  SNR.  As  in  the  case  for  the  N- 


correlator  receiver,  this  receiver  does  not  have  an  optimum 
threshold  that  is  obtained  from  the  theoretical  development. 
However,  a  procedure  similar  to  that  used  for  the  N-correlator 
receiver  has  been  used  to  obtain  the  optimum  threshold  setting. 
The  Pe  has  been  evaluated  as  a  function  of  SNR  except  that 
here,  plots  showing  Pg  as  a  function  of  y  have  not  been  included. 

Figures  4.13  through  4.16  have  been  obtained  through 
numerical  evaluation  of  Equations  (3.67)  and  (3.70).  Two 
curves  in  each  figure  have  been  plotted  to  show  the  effect 
of  the  parameter  m  on  the  receiver  performance.  Observe  that 
Equations  (3.67)  and  (3.69)  involve  a  parameter  a  which  can 
only  take  on  values  iT/N,  where  i  =  1,2, ...,N.  Observe  that 
i  =  N  corresponds  to  a  =  T  or  equivalently  a  =  0  since  v(t) 
is  T  periodic. 

Two  values  of  a  have  been  considered  for  each  of  the 
signals,  namely  the  sine  wave  and  square  wave,  in  order  to 
determine  its  effect  on  receiver  performance.  An  inspection 
of  Fig.  4.13  and  Fig.  4.14  displaying  the  performance  of  the 
receiver  for  v(t)  a  sine  wave,  or  inspection  of  Fig.  4.15  and 
Fig.  4.16  displaying  the  performance  of  the  receiver  for  v(t) 
a  square  wave,  it  can  be  observed  that  a  does  not  have  any 
significant  effect  on  Pg  for  either  case.  The  logic  support¬ 
ing  this  observation  on  receiver  performance  as  a  function  of 
a  follows  from  the  fact  that  a  is  just  the  estimate  of  the 
time  delay  of  the  signal,  which  is  obtained  by  choosing  that 


9: 


correlator  branch  that  produces  the  largest  output.  The 
actual  value  of  a  is  irrelevant  in  so  far  as  the  performance 
of  the  receiver  is  concerned.  The  receiver  simply  uses  the 
estimate  of  the  time  delay  to  produce  a  decision  based  on  a 
correlation  operation.  However,  one  of  the  factors  that  might 
improve  or  degrade  the  error  probability  is  the  accuracy  of 
the  estimate  of  the  time  delay  of  the  signal,  which  is 
itself  a  function  of  N  and  SNR. 
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PERFORMANCE  (SQUARE  WAVE) 


Figure  4.16  Performance  (Square  Wave),  (Alpha  = 


V.  CONCLUSIONS 


In  this  thesis,  the  analysis  and  design  of  noncoherent 
digital  receivers  that  optimally  (and  in  some  cases  subopti- 
mally)  detect  an  arbitrarily  shaped  periodic  waveform  with 
random  time  delay,  has  been  carried  out.  The  probability 
density  function (p.d.f . )  of  the  random  time  delay  \  is  assumed 
to  be  of  some  known  form.  A  uniform  p.d.f.  and  a  (parameter 
varying)  non-uniform  p.d.f.  have  been  considered  and  their 
effect  on  the  receiver  performance  has  been  studied. 

The  optimum  (in  minimum  probability  of  error,  P  ,  sense) 
receivers  have  been  designed  on  the  basis  of  likelihood  ratio 
test  under  low  SNR  assumptions.  The  design  of  suboptimum 
receivers  is  based  on  heuristic  approaches  that  intuitively 
yield  structures  or  algorithms  that  yield  good  performance. 

The  analysis  carried  out  in  this  thesis,  allowed  for  the 
possibility  of  unequal  prior  probabilities  of  occurrence  of 
the  hypotheses,  however  all  examples  and  computations  assumed 
equal  prior  probabilities  for  simplicity. 

The  low  SNR  assumption  allowed  certain  mathematical 
complexities  to  be  overcome.  Therefore  graphical  results 
involving  the  performance  (P  )  of  optimum  receivers  are  valid 
and  therefore  presented  only  for  SNR  values  less  than  or  equal 
to  -4  dB.  Operating  at  SNR  values  less  than  or  eaual  to  0  dB 
is  the  basic  requirement  of  covert  communication  and  may 
also  be  necessary  in  certain  radar  applications. 


The  optimum  receiver  operating  on  the  signals  that  have 
non- zero  DC  component,  detects  the  signal  with  higher  level 
of  error  probability  and  it  decreases  with  the  increase  of 
SNR  as  expected.  While  it  is  clear  that  no  receiver  could 
operate  with  such  an  undesirable  level  of  error  probability, 
the  graph  displaying  the  performance  (Pe)  does  give  an  indi¬ 
cation  of  the  performance  level  that  can  be  expected  at  low 
SNR's.  In  the  region  of  high  SNR,  it  is  expected  that  error 
probability  will  continue  to  decrease  to  reasonable  levels. 
However,  for  high  SNR's  this  receiver  is  no  longer  optimum. 

The  second  optimum  receiver  operating  on  arbitrary  sig¬ 
nals  does  also  yield  poor  performance.  But  here  again  we 
must  not  forget  that  the  receiver  is  operating  in  the  low 
SNR  region.  The  graphs  displaying  the  performance  of  this 
receiver  for  the  sine  wave  and  square  wave  give  an  indication 
of  the  performance  that  can  be  expected  at  low  SNR's.  Also, 
they  reveal  an  improvement  in  the  performance  with  an  in¬ 
crease  in  the  parameter  m.  It  is  concluded  that  for  high 
SNR's,  this  receiver  is  also  no  longer  optimum. 

One  of  the  receivers  designed  via  heuristic  means  was  a 
(suboptimum)  N-correlator  receiver.  Its  performance  was 
derived  and  tested  for  the  case  in  which  v(t)  was  a  periodic 
triangular  waveform.  While  this  receiver  is  unable  to  detect 
a  sine  wave,  it  is  capable  of  detecting  a  square  wave  or  some 
other  periodic  waveform.  The  other  receiver  designed  via 
heuristic  means,  namely  the  Estimator-Correlator  receiver, 
was  analyzed  in  terms  of  its  Pe  as  a  function  of  SNR. 


From  the  graphical  results  of  the  performance  of  the  N- 
correlator  receiver,  it  is  concluded  that  receiver  performance 
gets  better  as  N,  the  number  of  correlators,  gets  bigger. 

It  is  further  concluded  that  m,  the  parameter  that  determines 
the  uncertainty  about  the  time  delay  X ,  does  not  have  any 
effect  on  the  receiver  performance.  In  general,  it  was 
observed  that  except  for  the  case  of  the  N-correlator  receiver, 
the  parameter  m  had  a  significant  effect  on  the  performance 
of  the  receivers. 

While  a  number  of  practical  problems  have  been  addressed 
in  this  thesis,  related  problems  and  potential  studies  abound. 
For  example,  the  receiver  obtained  using  a  three-term  approxi¬ 
mation  in  the  exponential  of  Equation  (2.19)  needs  to  be 
analyzed  further  in  terms  of  performance  for  both  uniform  and 
non-uniform  p.d.f.  on  X.  Furthermore,  other  (suboptimum) 
receivers  need  to  be  proposed  and  analyzed,  such  as  simple 
energy  detection  receivers. 


APPENDIX  A 


DERIVATION  OF  MEAN  AND  VARIANCE  OF  R.V.  A 


The  p.d.f.  of  r.v.  A  (time  delay)  given  by  Equation  (2.1) 


is 


f  A  ( A  ;m) 


m  cos-^-A 

\  i0 757  x  i  lT/2l  ;■  0  i 


(A .  1) 


Substitution  of  Equation  (2.41)  into  Equation  (A.l)  yields 


f ^ ( A ;m) 


x  t  \  2tt\ 

Ip  (m)  cos  p-^-A 
tTq  (m) 


A  <  | T/2 | 


0  <  m  <  <» 


(A. 2) 


so  that  the  expected  value  or  mean  of  r.v.  A  can  be  expressed 
as 


T/2 

E  {  A  }  =  /  A  f  .  (  A  ;  m )  d  A 

-T/2 


TIn  (m)  4  Ip(ni) 

n  p=  —  oo  c 


■o 


T/2 

j  cos ( p2n A/T) dA 

-T/2 


=  0 


(A. 3) 
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since  the  product  of  X  (odd  function)  and  cos  (p2irX/T)  (even 
function)  is  an  odd  function  and  integral  of  an  odd  function 
produces  zero. 

The  variance  of  X  is  given  by 


T /2 

E(X}2  =  /  X2f.(X;m)  dX 

-T/2  A 


T/2 


TIq  (m) 


.  2  mcos2TiX/T  , . 
X  e  dX 


-T/2 


(A. 3) 


since  we  have  a  zero  mean  random  variable. 

The  variance  is  clearly  a  function  of  m.  However,  the 
integral  of  Equation  (A. 3)  is  not  readily  evaluated  in  closed 
form.  The  computer  has  been  used  to  evaluate 


Var{X}  1  (  2  _mcosTTx 

T2  ~  4l0(m)  Q>  X  e  dX 


(A. 4) 


2 

and  it  has  been  clearly  demonstrated  that  Var{X}/T  becomes 
smaller  as  m  grows. 
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INTERPRETATION  OF  A  SMALL  EXPONENTIAL 
IN  EQUATION  (2.22)  AND  EQUATION  (2760) 


In  Equation  (2.22)  and  Equation  (2.60)  it  has  been  assumed 

that 

T 

o  0  A 

—  /  r(t)v(t-A)dt  *  e  <<  1  (B.l) 

o  0 

This  appendix  is  devoted  to  analyzing  under  what  conditions  is 
this  assumption  valid.  Observe  that  e  is  a  r.v.  so  that 
Equation  (B.l)  must  be  analyzed  in  a  statistical  sense.  Recall 
that 


E{ e | HQ }  =  0  (B. 2) 

since  the  noise  has  been  assumed  to  be  zero  mean  and 

T 

E{ £ | H. }  =  J-  I  v2(t-X)  =  E  (B .  3 ) 

o  0  No 

where  E  is  the  signal  energy  for  0  £  t  £  Tg . 

Thus,  provided  E/Nq  <<  1,  £  under  either  hypothesis  on 
the  average  takes  the  value  zero  or  a  value  close  to  it.  The 
excursions  of  £  away  from  the  mean  can  be  obtained  by  evaluat¬ 
ing  the  conditional  variances.  Since 


E{e  |Hq} 


and 

E{e2|H1> 


Thus 


=  E{[^-  /  n(t;  .'(t-A)dt]^} 

No  QJ 


4  //  E{n(t)n(t)  }v(t-A) v(x-A)  dt  dx 

N  0 
o 


/  v  (t-A)dt 

0 


( 


=  E{  [ 


Tn 

■5  ■)  0  o 

v  (t-A}dt  ■+  rj=-  /  n  ( t)  v  (t-A )  dt]  z  } 
No  0 


E{(|^)2  +  2(4-  /  n(t)v(t-A)dt) 


N  0 
o 


+  (4-  /  n  (t)  v(t-A)  )  } 

wo  0 


, 2E> 2  2E 

'n  J  n 

o  o 


( 


Var{e|H  } 


2E 
N  _ 


Varie | H1 } 


Therefore,  provided  that  E/N0  <<  1,  e  remains  very  close 


to  its  average  value  which  is  either  zero  or  very  close  to 
zero.  Note  also  that  in  a  more  rigorous  sense,  we  can  bound 


P{e2  > 6  jHi>  < 


E{ e2 | Hi> 
6 


(B .  7 


If  we  set  6=1,  then 


p{£2  >1|HQ} 


< 


2E 

N 

o 


and 


P{e2  >1|H1} 


2e.  2 

V 


(B .  8 


(B .  9 


In  both  cases  if  2E/N  <<  1,  then  P{ e  > 1  or  e  < -1 }  is  also 

o  _  _ 

small.  Therefore,  a  small  E/NQ  guarantees  that  the  assumption 
of  a  small  exponential  in  Equations  (2.22)  and  (2.60)  is 
valid  with  a  high  degree  of  probability. 
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